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Abstract 

Abstract: The phase structure of self-avoiding polymerized membranes is 
studied by extensive Hybrid Monte Carlo simulations. Several folding tran- 
sitions from the flat to a collapsed state are found. Using a suitable order 
parameter and finite size scaling theory, these transitions are shown to be of 

first order. The phase diagram in the temperature-field plane is given. 
PACS numbers: 87.22.Bt - 05.70.Jk 
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Fluctuating membranes and surfaces are basic structural elements of biological systems 
and complex fluids. Recent theoretical work and experimental studies indicate, 
that these sheetlike macromolecules should have dramatically different properties than linear 
polymers. Polymerized membranes which contain a permanently cross-linked network of 
constituent molecules have a shear elasticity, giving them a large entropic bending rigidity. 
In the absence of self-avoidance, polymerized [gj and fluid |§ membranes adopt a crumpled 
random structure. Theoretical predictions by Flory mean-field approximation and Monte 
Carlo simulations [f| and renormalization group studies supported the existence of a 
high temperature crumpled phase for self-avoiding polymerized membranes also, suggesting 
a possible finite temperature crumpling transition in the presence of an explicit bending 
rigidity ||. However, more extensive computer simulations [3,PHTl|1 found no crumpling of 
self- avoiding tethered membranes in a good solvent. This prediction was confirmed recently 
by experimental studies of graphitic oxide ||. 

On the other hand, polymerized vesicles undergo a wrinkling transition Q], and upon 
addition of 10 vol % acetone, Spector and co-workers || found small compact objects, which 
appeared to be folded. A poor solvent leads to (short-ranged) attractive interactions, and a 
single membrane was found to be flat for high temperatures |9[], but in a collapsed state for 
sufficiently low temperatures 0. The transition between the flat and the collapsed states 
of the membrane proceeds through a sequence of folding transitions, which were first found 
by cooling of a single membrane from the flat phase |12j. Because no hysteresis was found, 



it was ruled out that the folded configurations are metastable states. However, this method 
does not give sufficient evidence of the order or even the existence of a transition. For 
instance, hysteresis can also be found at second order phase transitions of finite systems and 
the results for one system size may be misleading. In addition, the experimentally observed 
wrinkling transition is first order Q. 

On that account, we present a systematic finite size scaling analysis of the folding tran- 
sitions. The Hybrid Monte Carlo algorithm was used, which provides for simulations of 
the canonical ensemble with constant temperature. Moreover, ensemble averages are inde- 



pendent of the discretization step size St, i.e. systematic errors |13|,14]|. A suitable order 
parameter is defined and a phase diagram is presented. 

Besides the nearest neighbour interactions, the membranes were modeled similar to those 



of Abraham and Kardar ||12|| . The iV particles of the polymerized membrane form the sites 
of a hexagonal shaped triangular lattice. The bond potential between nearest neighbour 
particles is 

^ B = E(^o-^) 2 , (i) 

m 

with an equilibrium length 6 an d distances between particle i and its nearest neighbors 



In place of this harmonic potential, tethers were used by Abraham and Kardar |L2 
All particles interact through a shifted Lennard- Jones potential 



yu = i y\j t%) Vr c 12 rl) ' v ^ 

, r c < Tij 

with a cut-off at r c = 2.5 . The repulsive part of this interaction guarantees self-avoidance 
of the membrane. 

The folding of the membrane can be described by the eigenvalues Ai 2 >= A2 2 >= A3 2 
of the inertia tensor 

1 N 

T »,P = Tf Z ( r ^ r */3 - r « r /3) > ( 3 ) 
iV i=l 

where a,/3 G {x,?/,^}, and the sum runs over all particles of a given configuration; r a is 
the a component of the center of mass for a configuration. We can estimate the expected 
change of the membrane eigenvalues by those of an unfolded or folded disc with radius 1 and 
width d = 0.1 as shown in Table |. Taking into account all three eigenvalues, an unfolded, 
folded or twice folded configuration can be distinguished. A collapsed configuration would 
be indicated by approximately equal eigenvalues. 

Membranes with up to iV = 1141 particles were simulated at inverse temperatures in 
the range (3 G [0.2; 0.75]. Configurations at (5 — 0.20, 0.27 and 0.40 were found to be in an 
unfolded, folded and twice folded state by visualization (figure not shown). 



At the first folding transition, the eigenvalue stays constant approximately, while 
the second eigenvalue \\ decreases by a factor c m 0.2790. Therefore, we define an order 
parameter m by 

' 1 



(m((3)) 
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where \2{0o) is a reference value of an unfolded membrane at (3q far below the inverse criti- 
cal temperature (3 C . At a first order transition, the average order-parameter discontinuously 
jumps at (3 C [|i5|-|T7|. Contrary, at a continuous phase transition we expect a power law behav- 
ior of the average order parameter (m)((3) with an exponent /3', i.e. (m) cx {(3 — (3 C Y ,0>f3 c 
|18|-pT|. Indeed, the average order parameter (m) becomes very steep near f3 c with increas- 
ing number of particles N (data not shown). In principle, one can measure the increase of 
the slope and compare to the predictions of finite size scaling theory. Of course, the slope 
of (m) ((3) is proportional to the susceptibility x, which can be measured by the fluctuations 
of m also: 

X (f3) = L d f3 ((m 2 ) L - (mf L ) (5) 

The susceptibility xW) is shown in Figure [I] and the scaling of the maximum maxx(/3) in 
the inset of Figure [l|. 

At a first order transition, x{Pc S ) is expected to increase proportional to N = L d |T5|-[T7|. 
At a continuous phase transition, \(3 — /?c|~ 7± is predicted [pT^ -pl|. The L d increase of x{P) 
in Figure |I| gives evidence of the scaling at a first order transition. In addition, the width of 
the susceptibility peak should decrease as L~ d [|15| --|T7jl . Figure ^| shows a finite size scaling 
plot of the susceptibility data. Within errors, first order scaling behavior can be observed, 
at least above the inverse transition temperature. 

The transition temperature can be extrapolated by the position of the maximum in 
x{0) the minimum of the cumulant []T7| and the equal weight criterion |23| of the order 
parameter distribution. At least for x, the equal weight criterion predicts a shift of the 
effective transition temperature proportional to L~ 2d [22H24|, whereas a shift proportional 



to oc L~ l l u is expected for a continuous phase transition. Using quadratic terms in the 
regression, the extrapolations in Figure of the three observables agree within errors. The 
transition temperature of the first folding transition is found to be = 0.247(5). 

The increase of the susceptibility is caused by the characteristic double-peak structure of 
the order parameter distribution near f3 c , which is typical for a discontinuous phase transition 
T5|- p!8| , p2| -p5[] . Figure f| shows the expected double-peak distribution P{m) at the equal 



height transition temperature The development of a minimum in P{m) is confirmed 
by the method of Lee and Kosterlitz ||25|| . The measured AF in the inset of Figure || is 
proportional to the free energy difference at the equal height transition temperature and 
increases as (L d ) x , x~1.3. 

Further, P{m) can be described by the reduced cumulant Ul 

U L (3) = l--^h . (6) 

At a continuous phase transition, Ul(P) is expected to approach 2/3 for all (3. The data 
shows a minimum, which becomes more pronounced for large N, indicating a first order 
phase transition (figure not shown). 

The folding of the membrane must be visible in the attractive part of the potential 
energy, also. In fact, there is a jump in the potential energy and the related specific heat 
develops a peak, although very slowly. Besides the above defined order parameter, which is 
based on the geometry of the membrane, we can derive a different order parameter from the 
attractive part of the Lennard- Jones potential: 



m 



i{ge (rc - ry )(24 + 2 I)} . (7) 



Compared to m, the order parameter m has the advantage, that it is a local property. We 
performed Hybrid Monte Carlo simulations with an additional term h-fh in the Hamiltonian, 
where h is the conjugate field. Figure |5] shows the phase diagram in the (/?, fo)-plane of a 
membrane with 271 particles. The first order transition lines were computed by the multi- 
histogram method For the transition from one to two folds, an order parameter similar 



to m is used, which is based on Ai instead of A2. h = 1 is an upper limit for both transition 
lines, because of the vanishing attractive interactions at h — 1. 

The main result of this study is the first order nature of the folding transitions of self- 
avoiding polymerized membranes. This is in agreement with the observed folded structure 
of graphitic oxide |3[] and the first order wrinkling transition of partially polymerized vesicles 

I- 

Recently, another first order folding transition was found by DiFrancesco and Guitter 
27fl for regular triangular lattices embedded into two-dimensional space. However, these 



transitions combine flat and compact states directly, i.e. without several distinct foldings. 
Therefore, the folding transition described here is of a different kind. 

The unfolding of a singly folded membrane bears close resemblance to the unbinding 
transition of two distinct surfaces. Regarding a folded membrane of N particles, the fraction 
of particles near the crease decreases with 1/yN. Therefore, the nature of the folding 
transition is related to the nature of the underlying unbinding transition of two distinct 
membranes without the crease. The shape fluctuations of a single membrane of lateral size 
£|| are characterized by the typical fluctuation amplitude £j_ oc f;\\* . Polymerized membranes 
without lateral tension have a roughness exponent ( pa 0.6. The steric hindrance of two 
interacting membranes at separation I leads to an overall loss of entropy, which can be 
regarded as an effective fluctuation- induced repulsion, V rep oc 1// T with decay exponent 
r pa 3.3 for polymerized membranes. This repulsive interaction causes the unfolding of the 
membrane even in the presence of attractive van der Waals interactions. However, the crease 
of the folded membrane introduces an additional attractive interaction. This situation is 
similar to a membrane interaction which exhibits two minima at two different separations. 



Such an interaction implies a first- order unbinding transition |28| and may be an explanation 
of the first-order folding transition described in this work. 

Further investigations are necessary to determine the nature of the relation between 
folding and unbinding of polymerized membranes. 
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TABLES 



configuration 


Ai 2 


A 2 2 


A 3 2 


disc 


0.2543 


0.2473 


0.00084 


disc, folded 


0.2524 


0.0690 


0.00332 


disc, folded twice 


0.0929 


0.0487 


0.01330 


disc 


1.0000 


1.0000 


1.00 


disc, folded 


0.9925 


0.2790 


3.95 


disc, folded twice 


0.3653 


0.1969 


15.83 



TABLE I. Eigenvalues of the moment of inertia tensor of a disc with radius 1 and width d = 0.1. 
In the lower part, all numbers are relative to the values of the unfolded disc. 
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FIGURES 
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inverse temperature p 



FIG. 1. Susceptibility x(P) f° r membranes with N = 127,271,547,817 and 1141 particles. The 
solid lines were computed by the multi-histogram method. The scaling of the maximum maxx is 
shown in the inset. 



0.08 




iP-p;iL d 



FIG. 2. Finite size scaling plot of the susceptibility x(P) f° r a & TS ^ order transition using the 
value (3™ = 0.247(5) from Figure |. 
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0.000 0.001 0.002 0.003 0.004 

1 /L d 

FIG. 3. Critical temperature (3^° = 0.247(5) determined by the position of the maximal suscep- 
tibility x, the minimum of the cumulant Ul and equal-weight of the order parameter distribution 
for the four largest system sizes. 
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order parameter m 

FIG. 4. Probability distribution P(m, j3 ew ) of the order parameter m at the equal-weight transi- 
tion temperature (3 ew . The inset shows the free energy difference AF at the equal-height transition 
temperature (3 e h, which increases oc (L d ) x , xRil.3. 
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FIG. 5. Phase diagram of a membrane with 271 particles in the 0, h plane. The symbols denote 
the result of the simulations, the solid lines are the first order transition lines computed by the 
multi-histogram method. 
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